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Abstract 

We apply genus statistics to simulated CMB polarization maps, constructed from 
secondary-ionization cosmological models in experimental situations comparable to 
those of forthcoming space experiments. We find that both the cosmic baryon den- 
sity and the spectral index of density perturbations are strongly anticorrelated to 
the reionization redshift. Using the Fisher matrix approach we show that the accu- 
racies in determining the spectral index and the optical depth to the reionization 
epoch are better in the case of genus statistics than in standard power spectrum 
statistics. 
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1 Introduction 



The cosmic microwave background (CMB) offers one of the best probes of 
the early universe, enabling us to test the consistency of structure formation 
models. The detection of the CMB anisotropies at large scales by the COBE 
satellite (Smoot et al., 1992) as well as a number of new detections at inter- 
mediate and small scales [see e.g., White et al. (1994), Scott et al. (1995) and 
Bond (1996)] provided precious information to constrain several cosmological 
parameters. In recent years there has been an increasing interest for the CMB 
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polarization, leading to substantial progress regarding the structure and solu- 
tions of the transport equations (Ma & Bertschinger, 1995; Hu et al., 1998), 
numerical algorithms (Seljak & Zaldarriaga, 1996; Zaldarriaga et al., 1998), 
and statistical descriptors (Melchiorri & Vittorio, 1996; Seljak, 1997; Ng & 
Liu, 1998; Naselsky & Novikov, 1998). 

As CMB polarization is produced when the anisotropic radiation possessing 
non-zero quadrupole moment is scattered by free electrons via Thomson scat- 
tering (Rees, 1968), its magnitude, spatial distribution and topological proper- 
ties can provide information that is complementary to that obtained from the 
anisotropy alone. Unlike the temperature fluctuations which may evolve be- 
tween the last scattering and today, the CMB polarization probes the epoch of 
the last scattering directly. Also, different sources of temperature anisotropies 
(scalar and tensorial) give different signatures in the polarization power spec- 
tra. Therefore polarization is a powerful tool for reconstructing the sources 
of anisotropy. A full description of polarized radiation requires the introduc- 
tion of three additional power spectra, defining the "electric" and "magnetic" 
components of the polarization vector pattern and the cross correlation with 
the temperature anisotropies (Seljak & Zaldarriaga, 1996; Efstathiou & Bond 
1998). Including this additional information one can distinguish among phys- 
ical processes generating the temperature power spectrum, and hence better 
constrain cosmological models (Zaldarriaga, Spergel & Seljak, 1997). 

The current measurements of CMB temperature anisotropy already permit 
to place significant constraints on certain cosmological parameters. These pa- 
rameters include the amplitudes and the spectral indices of scalar and tensor 
perturbations (Knox & Turner, 1994; Crittenden, Coulson & Turok, 1995; 
Knox, 1995; Lidsey et al., 1997; Souradeep et al., 1998; Copeland et al., 1998), 
the various components of the mass density of the Universe and the Hubble 
constant (Jungman et al., 1996b; Lineweaver et al., 1997, 1998; Bond, Efs- 
tahiou & Tegmark, 1997; Zaldarriaga, Spergel & Seljak, 1997; Bond, Jaffe & 
Knox, 1998; Zaldarriaga, 1998; Hancock et al., 1998; Webster et al., 1998; 
Bartlett et al., 1998a, 1998b). The role of the CMB polarization in the deter- 
mination of the cosmological parameters has been discussed in several works 
(Zaldarriaga, 1997,1998; Seljak, 1997; Kamionkowski et al. 1997a; Efstathiou 
& Bond, 1998). Their results show that, although polarization can improve the 
accuracy of many cosmological parameters only by a modest amount, there are 
two important exceptions, regarding the discrimination of scalar and tensor 
modes and the determination of the reionization redshift of the intergalactic 
medium. Some degeneracies among cosmological parameters which are ex- 
pected to arise in future measurements of CMB anisotropies will be broken 
by the joint exploitation of polarization and temperature anisotropy (Kinney, 
1998), and others by the simultaneous use of some more astronomical data 
(Efstathiou & Bond, 1998). 
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The main limitation for the cosmological use of polarization is that it is pre- 
dicted to be small. The best experimental upper limit now available, being 
P ^ 16 /iK at a scale of 1.4° (Netterfield et al., 1995), is not very stringent 
from a theoretical point of view. However, a number of experiments which 
are being planned both at ground [POLAR (Keating et al., 1998), MITO (de 
Petris, 1998)] and from space [MAP (Bennet et al., 1996), SPOrt (Cortiglioni 
et al., 1997), PLANCK (Bersanelli et al., 1996)] are expected to improve the 
situation quite substantially. In view of the best exploitation of such forth- 
coming experiments (mainly, of MAP and PLANCK) detailed work has been 
performed on statistical descriptors of polarization. However, previous full-sky 
statistical studies of the polarization are generally based on the power spec- 
trum estimators and cross- correlations in Fourier (Seljak & Zaldarriaga, 1997; 
Kamionkowski et al., 1997a; Zaldarriaga & Seljak, 1997) and real space (Ng 
& Liu, 1998). A larger variety of methods have been proposed for studies of 
temperature maps: These include the analysis of 2-point correlation functions 
(Hinshaw et al., 1996), power spectra in terms of spherical harmonics mod- 
ified for Galactic cut (Wright et al., 1994) and in terms of Karhunen-Loeve 
eigenmodes (Bunn & Sugiyama, 1995), eventually including data compression 
(Tegmark et al., 1997), and topological methods (Torres et al., 1995). As to 
the latter methods, in the case of polarization fields, only general topological 
properties have been treated in the recent works of Naselsky & Novikov (1998) 
and Dolgov et al. (1998). 

In this paper we analyze the topological structure of the CMB polarized field 
employing genus statistics in order to study the degeneracies among param- 
eters of cosmological models. Our choice is motivated by the fact that genus 
is a locally invariant statistical estimator (Bond & Efstathiou, 1987; Gott et 
al., 1990; Torres et al., 1995; Schmalzing & Buchert, 1997), in the sense that 
an incomplete and non-uniform sky coverage leaves this quantity unchanged. 
In order to consider a realistic experimental situation, we perform simulations 
adopting the large but incomplete sky coverage consistent with the environ- 
ment on board the International Space Station Alpha (ISSA) and actually 
planned for SPOrt. We also adopt the angular resolution FWHM = 7° com- 
peting to POLAR and SPOrt, and to previous topological analyses of COBE's 
anisotropy field (Fabbri & Torres, 1996). 

The above angular scale is very interesting for cosmological models with a 
secondary ionization (Seljak & Zaldarriaga, 1997; Kamionkowski et al., 1997b; 
Zaldarriaga & Seljak, 1997), which predict the existence of a broad peak in 
the polarization power spectrum at low order multipoles (/ ^ 30) that is 
not present in the anisotropy power spectrum (Zaldarriaga, 1997). In flat- 
space models, the peak position scales as I oc rJ/^/if^Xe)" 1 / 3 , where r ri is the 
reionization optical depth, the density parameter of the baryonic matter, h 
the reduced Hubble constant and x e the ionization fraction. In general, large- 
scale polarization is enhanced in reionization models (Bond & Efstathiou, 
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1987; Zaldarriaga & Harari, 1995; Ng & Ng, 1995), the main peak amplitude 
being roughly proportional to r ri for r ri ^ 1 (and also dependent on the 
spectral index of the primordial power spectrum) (Zaldarriaga, 1997). For 
many of our model simulations we assume a rather strong reheating, T ri ~ 1, 
which is marginally consistent with experimental limits (de Bernardis et al., 
1997). In such a case we show that quite significant results can be obtained 
from genus analysis in experiments with a pixel sensitivity of 1 /iK, which is 
quite consistent with MAP's expected noise considering a 7° beam averaging. 
For weaker reionizations, which are more likely on experimental (de Bernardis 
et al., 1997; Fabbri, 1998) and theoretical (Haiman & Loeb, 1997a, 1997b) 
grounds, similar results can be obtained just scaling the noise level; as a matter 
of facts, the genus does not depend on the absolute amplitudes of signal and 
noise, but only on their ratio. Thus for instance, similar results for models 
with r ri ~ 0.1 can be obtained at the sensitivity level expected for PLANCK. 

The polarization genus test of cosmological models is affected by its own 
parameter degeneracies. Here we consider standard CDM models, assuming 
scalar modes (with primordial spectral slope n s ) and adiabatic initial condi- 
tions. Spanning the four-dimensional parameter space s 4 = (z ri , n s , Q b , x e ), we 
find that n s , Q b and x e are strongly anticorrelated to z ri . The high-likelihood 
regions in parameter space are found to be quite narrow from our simulations. 
Fixing z r i = 100 for both input and target models, the other parameters 
could be determined to a high accuracy for a 1 fiK detector noise; for in- 
stance, 5n s ~ ±0.05 and 5il b ~ ±0.008 at 68% CL. (All errors here and 
henceforth obviously include cosmic variance.) We conclude that a topological 
analysis of polarization can contribute to obtain a very accurate determination 
of parameters when combined with some other test. 

The most common way to evaluate cosmological parameters is to use the power 
spectrum estimates in connection with the standard likelihood function (see 
e.g. Bartlett et al., 1998b; Efstathiou & Bond, 1998; Zaldarriaga 1998). It is 
thereby worth making a detailed comparison of the advantages and disavan- 
tages of the genus technique versus the power spectrum technique. To this 
purpose, we computed formal errors by means of the Fisher information ma- 
trix method (e.g., Bond et al., 1997; Tegmark et al., 1998). This method allows 
us to quickly compare the relative accuracies of different techniques, although 
it cannot be used for a very accurate determination of likelihood contours in 
parameter space. We performed calculations for the experimental 

configurations pertaining to MAP, SPOrt and PLANCK, and found that the 
genus analysis is more efficient than the power spectrum analysis for determin- 
ing r r i and n s , while the opposite is generally true for Q b (with the exception 
of the SPOrt configuration). We also found that combining anisotropy and 
polarization data (for instance, performing the genus analysis on both maps) 
typically reduces the errors by a factor ~ 3. 
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The plan of the paper is the following. In the next Section we present the for- 
malism and method used for Monte-Carlo simulations of the Stokes parameter 
Q and U . The analysis of the polarization maps in terms of genus statistics 
is given in Section 3, and Section 4 presents the main results concerning the 
confidence regions of cosmological parameters obtained from simulations. In 
Section 5 we compare the accuracies in the estimates of the cosmological pa- 
rameters within the Fisher information matrix approach. 



2 Polarization maps 

The all-sky polarization maps were obtained following the formalism described 
in Zaldarriaga & Seljak (1997), by expanding the Stokes parameters Q and U 
in spin-weighted spherical harmonics with random amplitudes [cfr. also Sazhin 
& Benitez (1995)]. 



2.1 Stokes parameters 

The CMB polarization field can be described by a 2 x 2 temperature pertur- 
bation tensor T^-, in terms of which the Stokes parameters Q and U and the 
temperature anisotropy T are given by (Kosowsky, 1996; Bondi & Efstathiou, 



1987) Q = (T n - T 22 )/4, U = (T 12 )/2, T = (T n + T 22 )/4. (Circular polariza- 



tion is not necessary because it cannot be generated by Thomson scattering.) 
The combinations Q ±iU are quantities of spin ±2 and for a given direction 
in the sky ft can be expanded in spin- weighted spherical harmonics ± 2 Y/ m , 



(Q±iU)(n) = Xl a ±2,/m ±2Yim(n), 



(1) 



l,m 



while the temperature anisotropy is 



T(n) 



(2) 



l,m 



The linear combinations 



0-E,lm — 



{a-2,lm + a -2,/m)/2, 



(3) 



(4) 
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can be directly related with power spectra Cxi, which in the Gaussian theory- 
are rotational invariant quantities: 

< a *X,l'm'> a XM >= S a' S mm' C Xl, (5) 

< a* Tl , m ,,a EM >= S li: S mm 'Cci, (6) 

Here X stands for T,E or B denoting the temperature, electric-parity and 
magnetic-parity polarization modes respectively, and Cci describes the E-T 
cross correlation. (The B-T and B-E cross correlations vanish because of the 
opposite parities.) From Eq. (1) we obtain Q(n) and U{n) (Zaldarriaga, 1997) 

Qip) = ~ Yl\. a EM X hirn{n) + ia B M X Vm(n)], (7) 

Im 

U(n) = - Y^Wbm^M^) ~ ^£,; m X 2)/m (n)], (8) 

Im 

with 

X 1M (n) = ( 2 Y lm +_ 2 Y lm )/2 = ^/(2l + l)/4irF 1M (0)e tm t (9) 
X 2>lm {n) = ( 2 Y lm -_ 2 Y lm )/2 = v /(2/+l)/47rF 2iim (^)e^, (10) 

where the functions Fi i m and F 2> i rn can be calculated in terms of Legendre 
polynomials (Zaldarriaga, 1997; Kamionkowski et al., 1997a). The conditions 

X l,lm = X l,l~m, X 2 ,lm = - X 2,l-m, 

(lE,lm = a *E,l-m> a B,lm = a *B,l-mi (H) 

make Q and U real. 

For the computation of Q and U we need sets of random realizations of ax,i m 
consistent with assigned power spectra Cxi an d with the Gaussian distribu- 
tions accounting for cosmic variance. 

2.2 Monte-Carlo simulations of the Stokes parameters 

A given cosmological model only provides the power spectra. We considered 
standard CDM models without and with reionization (denoted by sCDM and 
srCDM, respectively), assuming only scalar modes with primordial spectral 
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slope n s and adiabatic initial conditions. We spanned the four-dimensional 
parameter space 

s 4 = {z ri ,n s ,Q b ,x e ) (12) 



fixing h = 0.5 and Y p = 0.24. As input models we assumed either sCDM 
or srCDM with z ri = 100, but for the Monte-Carlo simulated grids of target 
CDM models we investigated a fully 4-dimensional volume of S4. The four grid 
steps were 20, 0.05, 0.01 and 0.1, respectively, and interpolation was further 
performed for computation of the confidence regions (see next Section). 

All the relevant power spectra Ctj, Ce,i and Cc,i were obtained using the 
CMBFAST code developed by Seljak and Zaldarriaga (1996). 

For each realization the coefficients aE,i m were obtained following the proce- 
dure given in Zaldarriaga & Seljak (1997): For each multipole / we diagonalize 



the correlation matrix: Mi 



\Cci C 



, then we multiply the square roots 



El 



of the eigenvalues of Mi by a pair of Gaussian random numbers and rotate 
back to the original frame. Following this procedure we obtained random re- 
alizations of aE,i m satisfying the correct correlation properties. 

We constructed simulated maps of the Stokes parameters adopting a pixeliza- 
tion scheme of the "igloo" type (Crittenden & Turok, 1998). The sky region 
seen by a 7° beam of an experiment on ISSA (—51°. 6 < 5 < 51°. 6) was di- 
vided into rows with edges of constant latitude; each row was cut by constant- 
longitude lines, the angular distance between two neighbor pixels of constant 
latitude being Act = FWHM/ cos 5. Although the pixels obtained have un- 
equal trapezoidal shapes (becoming nearly rectangular only close to the Galac- 
tic plane), we have the following advantage, that the pixel edges defined by 
the spherical-coordinate frame allow a fast integration of the spin-weighted 
spherical harmonics. For each pixel we calculate Q(9i,fa) and U(9i,fa) (with 
6i = 7r/2 — 5 and fa = a), taking l max = 30. Equations (7) and (8) with the 
conditions (11), including the finite beamwidth and adding the detector noise 
contribution, become 



30 1 

Q(0i, fa) = - E E [Bi(a E ,im X i,im + a EM XliJ] + Ni, 

1=2 m=0 
30 I 

U(6i, fa) = -*EE [Bi{a EM X 2M - a* EM X* 2M )} + N h (13) 

1=2 m=0 

where Bi = exp [— a 2 l{l + l)/2] is the Gaussian window function of the beam, 
with a = 0.425 x FWHM and FWHM = 7°, and N { is a random realization 
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Fig. 1. The integrated genus per unit area for some reionization models. From top 
to bottom z r i = 200, 150, 100, 50, and 0. For all of the curves we take = 0.05, 
h = 0.5, Y p = 0.24, x e = 1 and n s = 1. 

Fig. 2. The \ 2 distributions obtained taking as a target model srCDM, and as input 
models srCDM (continuous line) and sCDM (dashed line). 

of the noise per pixel. The rms noise depends on the instrument sensitivity 
and the observing time for each pixel. The average value of the rms detector 
noise per pixel used in the simulations was 1 fiK. 



3 Genus analysis 

The specific signature of the CMB polarization field P = y/Q' 2 + U' 2 obtained 
for different underlying cosmological models was analyzed using the Euler 
characteristic of the field, equivalent to genus per unit area, under the as- 
sumption of random Gaussian primordial density perturbations. We determine 
the integrated genus per unit area G(p) above some threshold p (Naselsky & 
Novikov, 1998; Dolgov et al, 1998) as 

G{p) = N max + N mm - N sad , ( 14) 

where N max , N min and N sad are the number densities of maxima, minima and 
saddle points of P, respectively, above the threshold. The expectation values 
in terms of the curvature we have: 

<G(p)>= i© 2(p2 - 1)£ ~*' (15) 

where a and o\ are the spectral parameters (Bond & Efstathiou, 1987b) 

< Q 2 >=< U 2 >-- 

< QiQj >=< UiUj > 

and we set Qi = dQ/dxi and Ui = dU/dxi. For each random realization of 
the Stokes parameters we calculated the genus distribution for 30 values of p. 
Figure 1 

presents the integrated genus distributions < G(p) > obtained averaging over 
a set of 400 realizations of the Stokes parameters for different underlying 
reionization models. For such models we take il b = 0.05, h = 0.5, Y p = 0.24, 



0- 



Oij 2 , 



(16) 
(17) 
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Fig. 3. The confidence contour at 1 a level (continuous lines) in the z r i-n s plane for 
target reionized CDM models in four parameter space, when the input model was 
srCDM. The dashed line gives the best-fit curve. 

Fig. 4. The confidence contour at 1 a level (continuous lines) in the z r i-Qb plane for 
target reionized CDM models in four parameter space, when the input model was 
srCDM. The dashed line gives the best-fit curve. 

x e — 1 and n s — 1. For each set of Monte-Carlo realizations obtained for a 
given target model we calculate the x 2 estimator defined as 

30 30 

X 2 = EE(< G t9 ( Pi ) >-< G in ( Pi ) >)A-- 1 

i=l j=l 

( < G*(p,-) > - < G in ( Pj ) >), (18) 

where < G tg {p) > and < G m {p) > are the ensemble- averaged integrated genus 
for Monte-Carlo realizations of the target and input model respectively, and A^- 
is the covariance matrix of the Monte-Carlo realizations of the target model: 



1 -^realiz 

A« = t; E (G k ( Pl )-<G( Pt )>) 

realiz k=l 

( G\p 3 )- < G( Pj ) >). (19) 



Figure 2 presents two x 2 distributions obtained taking srCDM as target model, 
and srCDM and sCDM as input models. 



4 The confidence regions 



The confidence regions in four- dimensional parameter space s 4 = (z ri , n s , Qb, x e ) 
was obtained as constant x 2 boundaries at Xmin + <^xl (Press et al. 1992), for 
v — 4 and Axl = 4.72 and 7.78 at l-a and 2-a level respectively. 

Figure 3 presents l-a confidence interval obtained for the spectral index of 
the scalar modes Tig clS db function of reionization redshift z r i. The contour is 
obtained taking as input model srCDM with s 4 = (100, 1,0.05, 1). We found 
that n s and z„ are highly anticorrelated. The statistical fit for our range of 
parameters gives 

n s -l = -(7.12 ± 0.58) x 10" 5 ^/ 2 (20) 



at 95 % CL. The reason for this result is that the polarization amplitude at 
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the peak of the power spectrum is roughly proportional to the reionization 
optical depth 

r ri « 3.8 x KT 2 (x e Q b h) fr 1/2 (l + z ri ) 3/2 , (21) 

and also increases with the spectral index n s (i.e., enhancing small scale per- 
turbations). Thus the effect of decreasing r ri is compensated by changing the 
spectral index n s by an amount proportional to r ri oc zf( 2 . 

It is worth noticing that genus does not depend on the absolute amplitude 
of the polarization field, but only on its angular power spectrum. Thus for 
ideal, zero-noise experiments we do not expect Eq. (20) to be valid any more. 
However, since the detector noise has a quite different spectrum from the cos- 
mological signal, the signal amplitude is important in practice. Equation (20) 
should be regarded as an analog of equations describing the n s — Q r ms-ps au- 
tocorrelation found from COBE-DMR anisotropy data, including the angular 
correlation function (Seljak & Bertschinger, 1993) and topological analyses 
(Fabbri & Torres, 1996). 

Figure 4 presents the confidence contours in Q b -z r i plane at 68 % CL, obtained 
in the four- dimensional parameter space for h = 0.5. Here again we find a clear 
anticorrelation described by 

tt b = (7.55 ± 0.36) x 10~ 2 - (2.33 ± 0.12) x lO" 5 ^ 372 (22) 

with errors at 95% CL, and an analogous anticorrelation is found for the couple 
x e -z r i. Equation (22), too, is qualitatively interpreted by means of Eq. (21). 
It should be noted however that r r j is not constant along the curve defined by 
(22). As a matter of facts, although r ri is the most important parameter, the 
polarization field does also depend on other parameters. 

If z r i is fixed, then the other parameters are determined to a high accuracy: 
We get 5n s ~ ±0.05 and 5Q b ~ ±0.008 at 95% CL. Setting z ri equal to the 
input value, Eqs. (20) and (22) give displacements An s ~ —0.07 and Af2 6 ~ 
0.002 with respect to the input values. These numbers can be interpreted as 
estimates of bias; however, since they are comparable to errors and derived 
from the analytic best-fit curves, they should probably be regarded as upper 
limits on bias. 

The above accuracies refer to models with a rather strong reheating, r ri ~ 1, 
which provide polarized signals of several fiK, greater than the assumed pixel 
sensitivity of 1 fiK. In an analysis of this kind, such a sensitivity should be 
warranted over a large sky coverage (about 80% of the full solid angle in our 
simulations, i.e. over ~ 600 pixels). This requirement is somewhat beyond 
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several of the currently planned experiments. However, MAP's expected sen- 
sitivity of about 20 fiK at a scale of 0.3° scales just to 1 fiK considering a 
7° beam averaging. For weaker reionizations, we should consider that results 
rather similar to those presented here can be obtained just scaling the noise 
level, since genus does not depend on the absolute amplitudes of signal and 
noise, but only on their ratio. An optical depth r r j = 0.1 is quite realistic in 
the light of ionizing source models (Haiman & Loeb, 1997a, 1997b) and an 
analysis of anisotropy data (de Bernardis et al., 1997); in that case the noise 
level that we request on a 7° scale is of order 0.1 fjK in order to get similar 
accuracies. This is consistent with PLANCK's planned sensitivity at 100 GHz. 



5 Accuracy on the estimates of cosmological parameters 

In this section we investigate how measurements of the CMB anisotropies 
alone and anisotropies plus polarization can constrain the relevant cosmologi- 
cal parameters. The errors on the cosmological parameter estimates have few 
dominant components: 

• A nearly exact or "geometrical" degeneracy, that leads to nearly identical 
power spectra (Bond, Eftathiou & Tegmark 1997; Efstathiou & Bond 1998) 
provided we have identical matter content, primordial power spectra and 
angular size distance to the the last scattering surface. 

• The cosmic variance, that results from comparing a theoretical statistical 
distribution of observables with a finite distribution represented by the data. 
Cutting out parts of the sky, such as the Galactic plane, increases the cosmic 
variance by a factor approximately inversely proportional to the fraction of 
the sky sampled. 

• The method used (x 2 , maximum likelihood, the Fisher information matrix 
approximation) can also bias estimates of the cosmological parameters. 

The standard way to estimate parameters is the maximization of the likelihood 
function. If the likelihood function C for a given particular data set D is a 
multivariate Gaussian in D, then: 

C (s|D) oc ^ exph^Co^D], 

where Cov is the covariance matrix that embodies the cosmological parameter 
data set s and the noise. If C can be expanded to quadratic order about its 
maximum, then the accuracy with which the parameters in a given cosmolog- 
ical model can be reconstructed from the data set D can be obtained using 
the Fisher information matrix F^, whose elements measure the width and the 
shape of the likelihood function around its maximum (Bond, Efstathiou & 
Tegmark, 1997; Efstathiou & Bond, 1998): 
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Fij = -triAA.l Ai = Cov-'CD) — . 
The minimum error that can be obtained on a parameter is then given by 

5s t = ^fJ\ (23) 

depending not only on the experimental parameters data set, but also on 
the target and input cosmological models and the number of cosmological 
parameters involved in the computation. The data set D used to construct 
the Fisher information matrix can be either a sky map, the power spectrum, 
the cross-correlation function, or some topological descriptor like genus. 

5. 1 Power spectrum statistics 

The most common way to compute cosmological parameters is to use the 
power spectrum estimates (see e.g. Bartlett et al. ,1998b; Efstathiou & Bond, 
1998; Zaldarriaga 1998). This method still has several difficulties that could 
lead to biased estimates of the parameters (Wandelt, Hivon & Gorski, 1998b; 
Bartlett et al., 1998a,b). They arise mainly because Galactic cuts in the map, 
non-uniform sky coverage, anisotropic noise and systematic effects (like those 
induced by the foregrounds subtraction) make C\ non-Gaussian correlated 
quantities. In fact, even in Gaussian theories the Cj, that represents the vari- 
ances of individual spherical harmonic coefficients, are x 2 distributed. Hence, 
the standard likelihood method is not strictly applicable. 

If the temperature anisotropy power spectrum alone is used, then the Fisher 
information matrix given by equation (23) can be written as 

^Ef^(C^. (24, 



If both anisotropy and polarization power spectra are used, the Fisher infor- 
mation matrix reads as (Zaldarriaga & Seljak 1997; Zaldarriaga 1997) 

l X,Y ub * Ub 



where X and Y stands for T, E, C and B power spectra and Cov" 1 is the 
inverse of the covariance matrix. For the purpose of this work we assume only 
scalar modes. Then the relevant covariance matrix elements in Eqs. (24) and 
(25) are: 
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^^^^^^ 
Cov(C 2 El ) = + \ )fsky (C El + w^Bf 2 ) 2 , 
Cov{C 2 cl ) = 2 [C 2 cl + (C Tl + w- 1 BT 2 ){C El + w^Bf 2 )}, 

Cov(C Tl Cci) = (2l+\)f sky Cci{CTl + w ~ 1jB '" 2 )' 

Cov(C El Cci) = , 2 1W C C i(Cei + Wp'Bf 2 ). (26) 

Here we set w = J2c w c with the sum performed over detector channels, w c = 
(o ~c, P ix6 *c,pix)~ 2 (Knox 1995) and wp = 2w; also, Bf = J2 c B 2 iW c /w accounts 
for the beam smearing and B 2 d = e~ l ^ l+1 ^ ls is the Gaussian beam profile, 
l s = V8 \n2(6 c )j^ hm and f s k y is the fraction of the sky used in the analysis. 

Table 1 lists the experimental parameters of the various experiments that 
we considered in our calculations. In the Table the parameter set for MAP 
(Bennet et al., 1996) refers to the current, updated configuration. Also, we 
label by SPOrt-like the parameter set used in Section 4 for the computation 
of the confidence regions, that is generally consistent with the parameter data 
set of SPOrt-ISS, except for a sensitivity level somewhat better than currently 
achieved (Cortiglioni et al., 1997). For each experiment we only consider the 
frequency channels where the cosmological signal is not expected to be masked 
by the Galactic foreground. 

We assume as target model standard reionized CDM model (srCDM) (fl — 1, 
h = 0.5, Qb = 0.05, n s = 1, r ri — 1, x e — 1) and compute the derivatives of 
the power spectra with respect to n s , r ri and Qb using the truncated Taylor 
series expansion (Efstathiou & Bond 1998) 

Q( Si ) = Q(s ) + (^] As t . (27) 



Here s is the parameter data set of the target model and s is the parameter 
data set of the input models that differs within 5% of s . 

Table 2 presents l-a errors on the estimates of the relevant cosmological pa- 
rameters obtained from the anisotropy alone and anisotropy plus polarization 
for the experimental parameters listed in Table 1 and few values of f s k y - For a 
correct comparison with previous work we should consider that the expected 
errors depend on the number of model free parameters (4 in our case) . Larger 
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errors obviously arise when as many as ~ 10 free parameters are fitted (Jung- 
man et al, 1996a, 1996b). 



5.2 Genus statistics 



The expectation value of genus depends on the coherence angle of the field 
9 C = 2 1 / 2 (T /o- 1 (Gott et al, 1990), which can be written as: 

c EiKi + mi + i)(QBf + w-iy 1 ) 



One should note that the shape of the genus curve is fixed by the Gaussian 
random-phase nature of the field and its amplitude depends only the power 
spectrum. For a Gaussian random field Q 2 C is Gaussian distributed and the 
standard likelihood method is fully applicable. 

We use the statistics of the coherence angle to compute the errors on the 
estimates of the cosmological parameters for the experimental data sets listed 
in Table 1. The Fisher information matrices can be obtained from Eqs. (24) 
and (25) making the following substitutions: 



dQ 


del 


dC xl 


dsi 


dC xl 


dsi 


dQ 


del 


dC Yl 


d Sj 


dC Yl 


d Sj 



Cov(C xl C Yl ) 



del 



dC xl dC Y i 



de 2 

' Cov{C xl C Yl ). 



(29) 



Although the coherence angle does not depends on the fraction of the sky 
involved in the analysis, its covariance does so through the covariance matrix 
of the Q. 

Table 3 lists la errors on the estimates of the cosmological parameters obtained 
using the coherence angle statistics for few values of f s k y - The inspection of 
the results listed in Table 2 and Table 3 shows that: 

• Introducing polarization improves the accuracy on cosmological parameters 
from both power spectrum statistics and genus statistics, and as expected, 
the accuracy on r ri is improved most, generally by a factor ~ 3. 

• The genus statistics is more sensitive to n s and r ri than the power spectrum 
technique, while in most cases Qb is better determined in the case of power 
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spectrum. 

• At the largest beamwidth, which cuts off high order harmonics, the relative 
performance of genus statistics improves significantly. 

The second of the above results may seem somewhat strange at first sight. 
However, it can be interpreted by observing the topological analyses weight 
higher order harmonics more strongly than other current techniques, so that 
it is no surprise that the high-likelihood elongated "hills" in parameter space 
usually have somewhat different slopes (Fabbri & Torres, 1996). Thus couples 
of mutually anticorrelated parameters may have larger and smaller errors, 
respectively, for geometrical reasons. The third result, too, can probably be 
interpreted in terms of the weight given to high order harmonics, which con- 
trasts the loss of efficiency arising from the largest beamwidth. 



6 Conclusion 

The genus technique obviously suffers from parameter degeneracy problems 
as any other technique, no matter how sophisticated, that can be applied to 
the analysis of CMB data. However, since the anticorrelation curves are not 
identical to those arising from other techniques, the joint utilization of at least 
two techniques can significantly reduce the confidence regions of cosmological 
parameters. In this connection, CMB polarization is useful in order to provide 
independent data sets. 

A topological analysis of polarization maps can contribute to obtain accurate 
determination of some cosmological parameters that determine the spectral 
amplitude, in particular n s and z ri . Moreover it seems especially appropriate 
for large-beam experiments which are affected by the suppression of high order 
harmonics. 
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Table 1. The experimental parameters. 



i/(GHz) 9 fwhm a pix /W- 6 w^/W- 



MAP 


60 


21' 


12.1 


5.4 


(Bennet et al., 1996) 


90 


12.6' 


25.5 


6.8 


Planck LFI 


70 


14' 


3.6 


0.215 


(Mandolesi et al., 1998) 


100 


10' 


4.3 


0.156 




100 


10.7' 


1.7 


0.028 


PLANCK-HFI 


150 


8' 


2.0 


0.022 


(Puget et al., 1998) 


220 


5.5' 


4.3 


0.047 


SPOrt-like 


60 


7° 


1.0 


14.15 


(Cortiglioni et al., 1997) 


90 


7° 


1.0 


14.15 
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Table 2. 1 — a errors on the estimates of the cosmological parameters obtained 
from the power spectrum statistics. 







Anisotropy 


Anisotropy & 


pUldl lZfiLlUIl 


J sky 




1. 


0.7 


0.5 


1. 


0.7 


yj.o 




An y i n 3 

Oll s a 1U 


7.21 


8.52 


9.26 


3.19 


3.81 




MAP 


At v If! 2 

u i fee a ±yj 


4.84 


5.68 


6.92 


2.30 


2.67 






AOl v ID 3 


1.15 


1.41 


1.62 


0.96 


1.16 


1 




An V 

on s x ru 


2.97 


3.55 


4.21 


1.02 


1.22 


1 /M 
1.44fc 


PFANPK TFT 


At x in 2 


2.41 


2.88 


3.4 


0.78 


0.94 


1 1 1 
1.11 




AOi v in 3 


U.oy 


U.4U 


U.00 


U.ZO 


U.oo 






5n s x 10 3 


2.79 


3.33 


3.93 


0.75 


0.89 


1.05 


PLANCK-HFI 


<5r rec x 10 2 


2.28 


2.73 


3.23 


0.58 


0.69 


0.82 




tfftfc x 10 3 


0.38 


0.45 


0.53 


0.15 


0.18 


0.21 




5n s 


0.42 


0.51 


0.61 


0.11 


0.14 


0.17 


SPOrt-like 


St 


2.12 


2.51 


2.97 


0.66 


0.80 


0.94 




5Q b x 10 2 


2.89 


3.43 


4.10 


1.22 


1.46 


1.73 
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Table 3. 1 — a errors on the estimates of the cosmological parameters from 
genus statistics. 







Anisotropy 


Anisotropy & polarization 


J sky 




1. 


0.7 


0.5 


1. 


0.7 


0.5 




Oll s A 1U 


3.74 


4.56 


5.27 


1.31 


1.59 


1.82 


MAP 


At v ID 2 
u i fee a ±yj 


2.11 


2.47 


2.95 


0.76 


0.89 


1.02 




AOl v ID 3 


2.99 


3.61 


4.27 


1.41 


1.69 


1.99 




An V i n 4 
OTl s X 1U 


8.12 


9.71 


11.49 


2.62 


3.13 


3.70 


PFANPK TFT 


At V in 3 

C / A -LU 


4.93 


5.89 


6.97 


1.59 


1.90 


2.25 




AOi v in 3 


U.oU 


l.Uo 


l.ZO 


U.oz 


U.oo 


U.4D 




5n s x 10 4 


6.88 


8.22 


9.73 


2.21 


2.64 


3.13 


PLANCK-HFI 


<5r rec x 10 3 


4.26 


5.09 


6.03 


1.37 


1.74 


1.95 




tfftfc x 10 3 


0.81 


0.96 


1.14 


0.29 


0.34 


0.41 




5n s 


0.1 


0.12 


0.15 


0.035 


0.042 


0.051 


SPOrt-like 


St 


0.3 


0.35 


0.42 


0.096 


0.11 


0.13 




5Q b x 10 2 


1.81 


2.18 


2.58 


0.66 


0.78 


0.93 
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